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Abstract 

In this paper, we will use a quantum operator which performs the 
inversion about the mean operation only on a subspace of the system 
(Partial Diffusion Operator) to propose a quantum search algorithm 
runs in 0(y/ N/M) for searching unstructured list of size N with M 
matches such that, 1 < M < N . We will show that the performance 
of the algorithm is more reliable than known quantum search algo- 
rithms especially for multiple matches within the search space. A 
performance comparison with Grover's algorithm will be provided. 

1 Introduction 

Quantum computers @ |H1 E] are probabilistic devices, which promise to do 
some types of computation more powerfully than classical computers (HI • 
Many quantum algorithms have been presented recently, for example, Shor 
[TTj presented a quantum algorithm for factorising a composite integer into 
its prime factors in polynomial time. Grover [TO] presented an algorithm for 
searching unstructured list of N items with quadratic speed-up over algo- 
rithms run on classical computers. 
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Grover's algorithm inspired many researchers, including this work, to 
try to analyze and/or generalize his algorithm IU OH El EH E3- Grover's 
algorithm perfomance is near to optimum for a single match within the search 
space, although the number of iterations required by the algorithm increases; 
i.e. the problem becomes harder, as the number of matches exceeds half the 
number of items in the search space ^3] which is undesired behaviour for 
a search algorithm since the problem is expected to be easier for multiple 
matches. 

In this paper, using a partial diffusion operation, we will show a quan- 
tum algorithm, which can find a match among multiple matches within the 
search space after one iteration with probability at least 90% if the number 
of matches is more than one-third of the search space. For fewer matches 
the algorithm runs in quadratic speed up similar to Grover's algorithm with 
more reliable behaviour, as we will see. 

The plan of the paper is as follows: Section 2 gives a short introduction to 
quantum computers. Section 3 introduces the search problem and Grover's 
algorithm performance. Section 4 and 5 introduce the proposed algorithm 
with analysis on its performance and behaviour. And we will end up with a 
conclusion in section 6. 

2 Quantum Computers 
2.1 Quantum Bits 

In classical computers, a bit is considered as the basic unit for information 
processing; a bit can carry one value at a time (either or 1). In quantum 
computers, the analogue of the bit is the quantum bit (qubit ^Hj), which 
has two possible states encoded as |0) and |1); where the notation | ) is 
called Dirac Notation and is considered as the standard notation of states 
in quantum mechanics [7j. For quantum computing purposes, the states |0) 
and |1) can be considered as the classical bit values and 1 respectively. 
An important difference between a classical bit and a qubit is that the qubit 
can exist in a linear superposition of both states (|0) and |1)) at the same 
time and this gives the hope that quantum computers can do computation 
simultaneously (Quantum Parallelism). If we consider a quantum register 
with n qubits all in superposition, then any operation applied on this register 
will be applied on the 2 n states representing the superposition simultaneously. 
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2.2 Quantum Measurements 



To read information from a quantum register (quantum system), we must 
apply a measurement on that register which will result in a projection of 
the states of the system to a subspace of the state space compatible with the 
values being measured. For example, consider a two-qubit system \<f>) defined 
as follows: 



|0) = a|00)+/?|01)+ 7 |10)+(5|ll), (1) 

where a, j3, 7, and 5 are complex numbers called the amplitudes of the system 
and satisfy \a\ 2 + \/3\ 2 + \'~f\ 2 + \5\ 2 = 1. The probability that the first qubit 
of |0) to be |0) is equal to (|a| 2 + |/3| 2 ). If for some reasons we need to have 
the value |0) in the first qubit after any measurement, we must try some how 
to increase its probability before applying the measurement. Note that, the 
new state after applying measurement must be re-normalized so the total 
probability is still 1. 



2.3 Quantum Gates 

In general, quantum algorithms can be understood as follows: Apply a series 
of transformations (gates) then apply the measurement to get the desired 
result with high probability. According to the laws of quantum mechanics 
and to keep the reversibility condition required in quantum computation, 
the evolution of the state of the quantum system of size n by time t is 
described by a matrix U of dimension 2" x 2 n 



m = u\i>), (2) 

where U satisfies the unitary condition: UW = I, where U* denotes the 
complex conjugate transpose of U and I is the identity matrix. For example, 
the X gate (NOT gate) is a single qubit gate (single input/output) similar 
in its effect to the classical NOT gate. It inverts the state |0) to the state 
|1) and visa versa. It's 2x2 unitary matrix takes this form, 



X 



1 

1 



(3) 



and its circuit takes the form shown in Fig.(£Q). Notice that, from now on 
we assume that a horizontal line used in any quantum circuit represents a 
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qubit and the flow of the circuit logic is from left to right. For circuits with 
multiple qubits, qubits will be arranged according to the notation used in 
the figure. 



(a\0)+P\l)) 



X 



(/3|0> + a|l» 



Figure 1: NOT gate quantum circuit. 



Another important example is the Hadamard gate (H gate) which has 
no classical equivalent; it produces a completely random output with equal 
probabilities of the output to be |0) or |1) on any measurements. It's 2x2 
unitary matrix takes this form, 

H = T2U -i 

and its circuit takes the form shown in Fig. ((21). 



(4) 



x) 



Figure 2: Hadamard gate quantum circuit, where x is any Boolean variable. 

Controlled operations are considered as the heart of quantum computing 
[2j, Controlled-f/ gate is the general case for any controlled gate with one or 
more control qubit (s) as shown in Fig.fjHla). It works as follows: If any of 
the control qubits |cj)'s (1 < i < n — 1) is set to 0, then the quantum gate 
U will not be applied on target qubit i.e. U is applied on \t) if and only 
if all |q)'s are set to 1. The states of the qubits after applying the gate will 
be transformed according to the following rule: 

|q) -> |cj) ; 1 < % < n - 1 . s 

\t) -> \t CU ) = U c ^- C ^ \t) {0) 

where CiC2...c n _i in the exponent of U means the AND-'mg of the qubits 

Ci, C2, C n -\. 



H 



73(1°) + 



|1» 
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|ci) - 

\C2) - 

| Cn-l)- 
\t) - 



u 



ci) 


|ci) 


c 2 ) 


|c 2 ) 


c n-l) 


C n -1 


tcu) 


It) 



-©- 



|ci) 
|c 2 ) 

I C n -l) 

Few) 



a.Controlled-U 



b.Controlled-NOT 



Figure 3: Controlled gates where the back circle • indicates the control 
qubits, and the symbol © in part (b.) indicates the target qubit. 

If U in the general Controlled-f/ gate is replaced with the X gate men- 
tioned above, the resulting gate is called Controlled-iVOT gate (shown in 
Fig.((5Jb)). It works as follows: It inverts the target qubit if and only if 
all the control qubits are set to 1. Thus the qubits of the system will be 
transformed according to the following rule: 



\ci) — > Iq) ; 1 < i < n — 1 

|*) -> \tcN) = \t © CiC 2 ...C n _i) 

where Cic 2 . . . c n _i is the AND-ing of the qubits c\, c 2 , . 
classical XOR operation. 



(6) 

. . , c n _i and © is the 



3 Search Problem 

Consider a list L of N items; L = {0, 1, N — 1}, and consider a function 
/ which maps the items in L to either or 1 according to some properties 
these items shall satisfy; i.e. / : L — > {0, 1}. The problem is to find any 
« 6 L such that f(i) = 1 assuming that such % must exist in the list. It was 
shown classically that we need approximately N/2 tests to get a result with 
probability at least one-half. Let M denotes the number of matches within 
the search space such that 1 < M < N and for simplicity and without loss 
of generality we can assume that N = 2 n . Grover's algorithm was shown to 

solve this problem j3j in O {^J N/Mj . In |TH], it was shown that the number 

of iterations will increase for M > N/2 which is un desired behaviour for a 
search algorithm. To over come this problem it was proposed in ^3] that 
the search space can be doubled so the number of matches always less than 
half the search space and iterate the algorithm tx / A^2N/ M times so the 
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Measure 



O [y/N/M 

Figure 4: Quantum circuit for the proposed algorithm. 



algorithm still runs in O U/NjMy But using this approach will double the 

cost of space/time requirement. In the following section we will present an 
algorithm that can find a solution for M > N/2 with probability at least 
92.6% after applying the algorithm once. 

4 The Algorithm 

4.1 Iterating the algorithm once 



For a list of size N = 2 n , the steps of the algorithm can be understood 
as follows as shown in Fig.Q: 

1- Register Preparation. Prepare a quantum register of n + 1 qubits all in 
state |0), where the extra qubit is used as a workspace for evaluating 
the oracle U f 



|ty o ) = |O) 0n ®|O). (7) 

Register Initialization. Apply Hadamard gate on each of the first n 
qubits in parallel, so they contain the 2 n states, where % is the integer 
representation of items in the list: 



|Wi) = [E^ ® J) \W,) = [ -1= X; \i) ) ® |0> . (8) 



i=0 
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3- Applying Oracle. Apply the oracle Uf to map the items in the list to 
either or 1 simultaneously and stores the result in the extra workspace 
qubit: 



N-l , N-l 



w *) = -mY, (i«'> ® i° © /«» = ^ E ® (9) 



2=0 1 = 

4- Partial Diffusion. In this step, we will define a new operator: Partial 
Diffusion Operator (Y) which works similar to the well known Diffusion 
Operator used in Grover's algorithm ^2] except that it performs the 
inversion about the mean operation only on a subspace of the system as 
follows: The diagonal representation of the partial diffusion operator 
Y when applied on n + 1 qubits system can take this form: 



Y = # 0n ®/(2|O) (0| - I)H® n ®I, (10) 

where the vector |0) used in Eqn. (jl(Jj) is a vector of lenght P = 2N = 
2 n+1 . Applying Y on a general system Ef^o* \k); where, \5k\ 2 = 1, 
can be understood as follows: Without loosing of generality, the general 
system can be re-written as, 

P-l N-l N-l 

E = E <*> M) ® + E Pi to) ® i 1 ))' ( n ) 

fc=0 j=0 j=0 

where {a, = Sk ■ k even} and = 5k : k odd}, then applying F on 
the system gives, 



Y[Eh \k)) = (H® n ® / (2 |0) (0| - /) ® /) E^fc |*) 

,fc=0 / fc=0 

= 2 (g> / |0> (0| H® n 7)^4 |/c) -^4 

fc=0 fc=0 
AT-1 iV-1 

= E (2 (a> - «,) (|j> ® |0» - E & (b'> ® 



(12) 
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Figure 5: Quantum circuit representing the Partial Diffusion Operator Y 
over n + 1 qubits. 



where (a) = Ylj=o a j i s the mean of the amplitudes of the subspace 
Yl'j^o a j ® i- e - a PPlyi n S ^e operator Y will perform the in- 
version about the mean only on the subspace; ^2^=0 a i (U) ® 1^)) an d 
will only change the sign of the amplitudes for the rest of the system; 
Yl!j=o Pj ® a chcuit implementation using elementary gates 
[2] is shown in Fig.©. 



The main idea of using the partial diffusion operator in searching is 
to apply the inversion about the mean operation only on the subspace 
of the system which includes all the states which represent the non- 
matches and half the number of the states which represent the matches 
while the other half will have the sign of their amplitudes inverted to 
the negative sign preparing them to be involved in the partial diffusion 
operation in the next iteration so the amplitudes of the matching states 
get amplified partially each iteration. The benefit of this is to keep half 
the number of the states which represent the matches as a stock each 
iteration to resist the de- amplification behaviour of the diffusion oper- 
ation when reaching the turning points as we will see when examining 
the performance of the algorithm. 

Let M be the number of matches, which makes the oracle Uf evaluate to 
TRUE (solutions); such that 1 < M < N; assume that indicates a 
sum over all % which are desired matches [M states), and £\ " indicates 
a sum over all % which are undesired items in the list. So, the system 
| W2) shown in Eqn.([5|) can be written as follows: 



S 



N-l 1 AT-1 

= -7= E " (i«'> ® + -tt? E ' (io ® i 1 ))- ( 13 ) 



1 = v 2=0 

Applying Y on IW2) will result in a new system described as follows: 

N-l N-l N-l 

m) = Bl e " (io ® io» + &i E ' (io ® io» + ci e ' (io ® 

i=0 i=0 i=0 

(14) 

where the mean used in the definition of partial diffusion operator is, 

, v fN-M\ 



and ai, 61 and ci used in Eqn. ()14j) are calculates as follows: 

a l = 2{a 1 ) — ; b x = 2 (an) ; ci = -=. (16) 

V -/V v iV 

Such that, 



(JV - M) a? + M6? + Mc? = 1. (17) 

Notice that, the states with amplitude 61 was with amplitude zero 
before applying Y. 

5- Measurement. If we measure the first n qubits after the first iteration 
(q = 1), we will get the desired solution with probability given as 
follows: 

i- Probability to find a match out of the M possible matches; 
taking into account that a solution \i) occurs twice as: (|z) (g> |0)) 
with amplitude 61 and ® |1)) with amplitude C\ as shown in 
Eqn.([14|). can be calculated as follows: 
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P.™ =M(6? + c?) 

= ^((%r) 2 +(^) 2 ) US) 
= 5(M)-8(#) 2 + 4(M) 3 . 

ii- Probability Pis to find undesired result out of the states can be 
calculated as follows: 

P^ = {N-M)a\. (19) 
Notice that, using Eqn.(|17|). 

P« + PW = 1. (20) 

4.1.1 Performance after Iterating the Algorithm Once 



n, where N = 2 n 


Max. prob. 


Min. prob. 


Avg. prob. 


2 


1.0 


0.8125 


0.875 


3 


1.0 


0.507812 


0.937500 


4 


1.0 


0.282227 


0.968750 


5 


1.0 


0.148560 


0.984375 


6 


1.0 


0.076187 


0.992187 



Table 1: First iteration performance with different size search space. 

Considering Eqn. (jl4j) and Eqn. (jl8)l we can see that the probability to find a 
solution varies according to the number of matches M in the superposition. 

From Tabled we can see that the maximum probability is always 1.0, 
and the minimum probability (worst case) decreases as the size of the list 
increases, which is expected for small M because the number of states will 
increase and the probability shall distribute over more states while the av- 
erage probability increases as the size of the list increases. It implies that 
the average performance of the first iteration of algorithm to find a solution 
increases as the size of the list increases. 
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To verify these results, taking into account that the oracle Uf is taken 
as a black box, we can define the average probability of success of the first 
iteration of algorithm; average(Ps^), as follows: 



N 

average(P^) = ^ £ n C m pP 

M=l 

= W ^2 M\(N-M)\ (Pi + Cl) ^1) 
N 

= whp (m-dkLm). ( im2 - 16MN + 8M2 ) 

where n Cm = m\(n~m)\ ls ^ ne num ber of possible cases for M matches. 

We can see that as the size of the list increases (N — > oo), average{P s ) 
shown in Eqn. (j21|) tends to 1. 

Classically, we can try to find a random guess of the item, which repre- 
sents the solution (one trial guess), we may succeed to find a solution with 
probability p^ dassical ) = M/N. The average probability of success can be 
calculated as follows: 

average{P^ dasstcal) ) = £ £ N classical) 



M=l 

N 

1 V- N.M 



N --- (22) 

2^ M\(N-M)\N 



E 

M=l 
1 

2" 



It means that we have an average probability one-half to find or not to 
find a solution by a single random guess even with the increase in size of the 
list. 

Similarly, Grover's algorithm has an average probabilty one-half after 
arbitrary number of iterations as we will see. It was shown in j3] that the 
probability of success of Grover's algorithm after qc iterations is given by: 



Pj 9G) = sin 2 ((2g G + 1)9), where, < 9 < - and sin 2 (fl) = — . (23) 

The average probability of success of Grover's algorithm after arbitrary 
number of iterations is as follows (Appendix A in [T5]): 
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average(P^) = ^J2 ^ sin 2 ((2g G + 1)6) = -. (24) 

M=l 

Comparing the performance of the first iteration of the proposed algo- 
rithm, first iteration of Grover's algorithm and the classical guess technique, 
Fig.© shows the probability of success of the three algorithms just men- 
tioned function of the ratio (M/N). 




0.2 0.4 0.6 0.8 1 



M/N 

Figure 6: A plot of the probability of success of the first the iteration of 
proposed algorithm P s , first iteration of Grover's algorithm Pi* '* and the 
classical guess p^ classica ^ as a function of the ratio (M/N). 

We can see from Fig.© that the probability of success of the first iteration 
is always above that of the classical guess technique. Grover's algorithm 
solves the case where M = N/4 with certainity and the proposed algorithm 
solves the case where M = N/2 with certainity. The probability of success 
of Grover's algorithm will start to go below one-half for M > N/2 while the 
probability of success of the proposed algorithm will stay more reliable with 
propabilty at least 92.6%. 
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4.2 Iterating the Algorithm 

If we consider iterating the algorithm, the iterating block will be applying the 
oracle Uf and the operator Y on the system in sequence as shown in Fig. (J3J) . 
To understand the effect of each iteration on the system, we will trace the 
state of the system during the first few iterations. Consider the system after 
first iteration shown in Eqn. (jl4j) . second iteration will modify the system as 
follows: 

Apply the oracle Uf will swap the amplitudes of the states which represent 
only the matches, i.e. states with amplitudes b\ will be with amplitudes c\ 
and states with amplitudes C\ will be with amplitudes b\ so the system can 
be described as, 

N-1 N-1 N-1 

m = «i£" go ® io» + d J2 ' go ® + & i E ' (io ® i 1 ))- ( 25 ) 

i=0 i=0 i=0 

Applying the operator Y will change the system as follows, 

N-1 N-1 N-1 

m = a 2 j2" (io ® io» + c 2 E ' (i«'> ® I 1 ))- ( 26 ) 

i=0 i=0 i=0 

where the mean used in the definition of partial diffusion operator is, 

(a 2 ) = ^({N-M)a 1 + Mc 1 ), (27) 
and d2, b 2 and C2 used in Eqn. (J26J) are calculated as follows: 

a 2 = 2 (a 2 ) - ai; 6 2 = 2(a 2 )-ci; c 2 = (28) 
and the probabilities of the system are, 



Pj 2) = M (bj + c|) = M (bl + 6?) . 

= M (a 2 2 ) = M(b 2 + c 2 f = M{b 2 - h) 2 . 
In the same fashion, third iteration will give the following system, 



(29) 



N-1 N-1 N-1 



U f \W 5 ) = I W 6 ) =a 2 J2" (10 ® |0» + c 2 E ' (10 ® l«» + &2 E ' (10 ® I 1 )) 

(30) 



i=0 «=0 i=0 



13 



N-l N-l N-l 



Y \W e ) = \W 7 ) = a 3 ^ " (l*> ® |0» + h £ ' (10 ® I )) + c 3 E ' (10 ® I 1 ))' 

(31) 



i=0 i=0 i=0 



where the mean used in the definition of partial diffusion operator is, 

(a 3 ) = ^((N-M)a 2 + Mc 2 ), (32) 
and a 3 , b 3 and c 3 used in Eqn. fJBTj) are calculated as follows: 

a 3 = 2 (a 3 ) - a 2 ; 6 3 = 2 (a 3 ) - c 2 ; c 3 = -b 2 , (33) 
and the probabilities of the system are, 

PP = M(bl + 4) = M(bl + bl). 

P$ = M (a|) = M (6 3 + c 3 ) 2 = M (63 - 6 2 ) 2 • 

In general, the system after q > 2 iterations can be described using the 
following recurrence relations, 

JV-l N-l N-l 

\W {q) ) = (10 ® |0» + ^£'(10® |0» + £ ' S |1», (35) 

i=0 i=0 i=0 

where the mean to be used in the definition of the partial diffusion operator 
is as follows: For simplicity, let y — 1 — ^ and s = then 

(«<?) = (ya ff -i + (1 - Z/K-i) , (36) 
and a q , b q and c 9 used in Eqn. (j35J) are calculated as follows: 

a q = 2 (a q ) - a 9 _i; a = s, ai = s (2y - 1) , (37) 

6q = 2 (a q ) - c,_i; 6 = s, 61 = 2sy, (38) 

c 9 = -feg-i; c = 0, ci = -s, (39) 

and the probabilities of the system are, 



14 



= M (b\ + c 2 ,) = M + bU) . (40) 

PW = M (a 2 .) = M (6, + c,) 2 = M (6, - 6 9 _ x ) 2 . (41) 
For q > 2, Eqn.tjHZj), Eqn. (j3~%j) and Eqn. (j3~9"j) could be re- written as follows: 

a q = 2ya q - 1 - a q _ 2 ] a = s; a 1 = s(2y-l), (42) 
b q = 2yb q _ 1 - b q _ 2 ; b = s; b x = 2sy, (43) 

c q = -b q -i] c = 0; Ci = -s, (44) 

Solving the above recurrence relations (In Appendix A), the closed forms 
are as follows (Proved in Appendix B): 

/ sin ((g+ 1)0) sinM)\ 
^ = S l sin ((9) '^inWj' (45) 

k = s i^m^) , (^) 



sin (0) 
/ sin (g^) 



(47) 



V sin (0) 

where y = cos (0) and < < |. The above closed forms can be expressed 
via the Chebyshev polynomials of the second kind U q (y) ^1] , which are 
defined as follows, 

sin ((g+ 1)0) 

= sin(6>) ' (48) 

This allows us to re-write the above closed form in terms of Chebyshev 
polynomials of the second kind as follows, 

a q = s (U q - £/,_!) , (49) 
b q = sU q , (50) 
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Cg = SU q -l, 

And the probabilities of the system, 



(51) 



P W = (i 



cos^^ + C/^). 



(52) 



P<£ = cos ((9) (Z7, 




(53) 



Such that 



+ p$ = M (fe 2 + c 2) + (AT - M) a 2 
= N(b 2 q + c 2 q )+2(N-M) c q b q 

= (sin 2 ((j + 1) 0) + sin 2 - 2 cos (0) sin ((j + 1) 0) sin 

= ( cos2 C#) sin2 ( ) ~ sin2 W cos2 (8) + sin2 (^)) 

= sMC 1 - sin2 (#)) sin2 (*) - sin2 (#) ( x - sin2 (*)) + sin2 Of )) 



4.2.1 Performance of Iterating the Algorithm 

Now, we have to calculate how many iterations, q, are required to find the 
matches with certainty or near certainty for different cases of 1 < M < N. To 
find a match with certainty on any measurement, then P^ must be as close 
as possible to certainty. To calculate the number of iterations, q, required to 
satisfy this condition, we need the following theorem. 

Theorem 4.1 Consider the following relation, 



where U q (y) is Chebyshev polynomials of the second kind, y = cos (6) and 
< 6 < f , then, 



sin 2 (6») 
sin 2 (6>) 
sin 2 (0) 



(54) 



(1 - cos + ££_!) = !, 



(55) 



71-6 



or 6 



7i 



q = 



26 



2 
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Proof From the definition of U q shown in Eqn. then Eqn. \55\) can take 
this form, 

( i-co SW) r^(fa±M + ^ l= i, 



sin 2 (0) sin 2 (0) 



or, 



sin 2 ((q + 1) 0) + sin s 



1 + cos (0) . 



Using simple trigonometric identities, the above relation may take the 
form, 

cos {2q6 + 20) + cos (2q0) + 2 cos (0) = 0. 
Using the addition formulas for cosine we get, 

2 cos (2q8) cos 2 (0) - 2 cos (9) sin (2q8) sin (9) + 2 cos (9) = 0, 
2 cos (9) (cos (2g0) cos (9) - sin (2g0) sin (0) - 1) = 0, 



cos (0) (cos(2g0 + 0) - 1) = 0. 

From the last equation we get, 

cos (0) = or cos (2g0 + 0) = cos (— 7r) , 

which gives the required conditions, 

7T it — 9 

9= 2 orq = ^r- 

Using the above result, and since the number of iterations must be integer, 
following the same fashion as shown in 0] , then the required number of 
iterations is, 



7T 



N 



2\[2 V M 



(56) 



where [ J is the floor operation. The algorithm runs in O (y/N/Af) with 
no contradiction with the prove of optimality shown in EI] . 
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5 Comparison with Grover's Algorithm 



First we will summarize the above results from both Grover's and the pro- 
posed algorithm before starting the comparison. The probability of success 
of Grover's algorithm as shown in j3] is as follows: 



where sin 2 (6) 



M . 

N ' 



Pj 9G > = sin 2 ((2g G + l)0), 
< < | and the required qc is, 



<1g 




(57) 



(58) 



For the proposed algorithm, the probability of success is as follows, 



p(D = (i _ cos (0)) 



/ sin 2 ((g + i; 
V sin 2 (9) 



+ 



sin 



sin 2 (0) 



(59) 



where cos (&) = ! — |f ; < # < ~, and the required q is, 



Q 




(60) 



Fig.(J7J) shows the probability of success as a function of the ratio M/N 
for both algorithms, after 2, 3, 4 and 5 iterations. It is clear from the graphs 
of the proposed algorithm that the probability will never return to zero once 
started and the minimum probability will increase as M increases because of 
the use of the partial diffusion operator which will resist the de- amplification 
when reaching the turning points as explained in the definition of the partial 
diffusion operator, i.e. the problem becomes easier for multiple matches, 
where for Grover's algorithm, the number of cases (points) to be solved with 
certainty is equal to the number of cases with zero-probability after arbitrary 
number of iterations. 

Another way to understand the behaviour of both algorithms is to plot 
the probability of success using the calculated number of iterations for each 
algorithm, Fig.© shows the probability of success as a function of the ratio 
M/N for both algorithms by inserting the calculated number of iterations 



qc and q shown in Eqn.(|3H|l and Eqn. (JoT)j) in Pg^ and P { s qGl respectively. We 
can see from the plot that the minimum probability that Grover's algorithm 
may reach is approx. 17.49 % when M/N = 0.61685 while for the proposed 



>(<?g) 
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Figure 8: Probability distribution using the appropriate number of iterations 
for both algorithms. 
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Figure 9: Probability distribution using the appropriate number of iterations 
for both algorithms for the hardest cases whereM / N < 1 x 10~ 3 . 
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algorithm, the minimum probability is 84.72% when M/N = 0.30842. Fig.© 
shows the same behaviour for both algorithms for small M and large N (hard 
cases where M/N < 1 x 10~ 3 ). It is interesting to notice that the behaviour 
for the proposed algorithm shown in Fig.fjHJ) is similar to the behaviour of 
the first iteration shown in Fig.© for M/N > 0.30842 which implies that if 
M/N > 0.30842 then the proposed algorithm runs in 0(1), i.e. the problem 
is easier for multiple matches. 



6 Conclusion 

In this paper, we presented a quantum algorithm for searching unstructured 

list of size N runs in O (^\/N/m\ where M is the number of matches within 

the list. The algorithm operations based on the Partial Diffusion Operator 
works similar to the Diffusion Operator used in Grover's algorithm JU] except 
that it performs the inversion about the mean only on a subspace of the 
system. Using this operator, we showed that the algorithm performs more 
reliable than Grover's algorithm in case of fewer number of matches (hard 
cases of the problem) and runs in 0(1) in case of multiple matches (easy 
cases of the problem). 
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Appendix A 

In this appendix, we will solve the following recurrence relations shown in 
Eqn. (J42j) and Eqn. (J43|) to get their close forms: First for a q which is defined 
as follows for q > 2, 

a q = 2ya q _i — a g - 2 ] q>2 such that ,a = s; a\ = s{2y — 1). (61) 

The characteristic equation, 

A 2 - 2y\ + 1 = 0, (62) 

and, 

\ 1 , 2 = y±iy / l-y 2 . (63) 
Let y = cos (9) such that < 9 < |, then, 

A 12 = cos (9) ± i sin (9) = e ±id . (64) 
So, the closed form will take the form, 

a q = Vl e iqe + v 2 e~ iqe , (65) 

where V\ and i; 2 are constants to be determined from the initial conditions 
as follows, 

s (e~ ie - 2y + 1) s (2y - 1 - e ie ) 

Vl = e -ie _ e e ' v * = ^Te _ e e • ( 66 ) 

Substituting in Eqn.(|65jl we get, 

a = -1— (2 cos (9) sin (q9) - sin ((q -1)9)- sin (q0)) . (67) 
sm (0) 
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We have from the identities of multiple-angle formulas, 

sin ((q + 1) 9) = 2 cos (9) sin (qO) - sin ((q -1)9). (68) 
So the closed form of a q may take the following form, 



/ sin ((g + 1)0) shH 
9 ~ V sin (0) sin (0) J • lWj 

Second, for 6 g which is defined as follows for q > 2, 

6 9 = 2yb q ^\ — b q _2', such that, b = s; bi = 2sy. (70) 

Since we are starting with the same recurrence relation used for a q but 
with different initial conditions, then the closed form for b q may take the 
form, 

b q = Vl e iqe + v 2 e~ iq \ (71) 

where V\ and v 2 are constants to be determined from the initial conditions 
as follows, 

se~ td — 2ys 2ys — se ld . . 

Substituting in Eqn. (|7T|) we get, 

\ = — (2 cos (9) sin (qO) - sin {(q - 1) 9)) . (73) 

So the closed form of b q may take the following form, 

/ sin 

~ sin (9) 

Appendix B 

In this appendix, we want to prove that, for q > 2, if the amplitudes a q , 
b q and c q used in Eqn. (|H5|) are defined according to the definition of Partial 
diffusion operator shown in Eqn. (fTU|) as follows, 

(a q ) = (ya,_i + (1 - y)c q -i) , (75) 



6g= / sin ((g + 1)6 (w) 
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a q = 2 (a q ) - a q -x\ a = s, a x = s (2y - 1) , 



(76) 



bg = 2 (a q ) - c ? _i; b = s, b x = 2sy, (77) 

c q = -bg-!] c = 0, ci = -s, (78) 

then their closed forms follows, where y = cos (0) and < 9 < |: 

Proof Substituting y = cos (6) in the definition and eliminating c q since it 
is sufficient to prove the closed forms for a q and b q we get, 

a q = (2 cos (9) - 1) a 9 _i - 2 (1 - cos (9)) b q . 2 , 



b q = 2 cos (9) a 9 _i — (1 — 2 cos (9)) b q -2, 
with initial conditions, 

a — s; ai = s(2cos(6 l ) — 1), 

b = s; &i = 2s cos (#) , 

Step 1: Prove for q = 2. 

For ci2, from definition and initial conditions, 
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a 2 = (2 cos (0) - 1) ai - 2 (1 - cos (0)) 6 

= (2 cos (0) - 1) (2 cos (5) - 1) s - 2 (1 - cos (0)) s 

= s ((2 cos (9) — l) 2 — 2 (1 — cos (0))) 

= s(4cos 2 (0) -2cos(0) - 1) 

= s (3 cos 2 (0) - sin 2 (0) - 2 cos (0)) 



For 62, from definition and initial conditions, 



b 2 = 2 cos (0) (2 cos (0) - 1) s - (1 - 2 cos (0)) s 
= s(4cos 2 (0) - 1) 
= s(3cos 2 (0) - sin 2 (0)) 



Step 2: Assume the relation is true for q = t — 1 and q — t, 



S V sin (0) sin (0) J ' a< S \ sin (0) sin (0) 



Step 3: Prove for g = £ + 1, 

For Of+i, from the definition and assumption, 

a m = (2 cos (0) - 1) at - 2 (1 - cos (0)) b t ^ 



= (2 cos (0) - 1) s {^^1 - Sg) - 2 (1 - cos (0)) s (0) 
= (2 cos (0) sin ((t + 1) 0) - sin {t9) - sin ((t + 1) 0)) 



For bt+i, from the definition and assumption, 









2(3 



b t+1 = 2 cos (9) a t -(l-2 cos (0)) 

= 2 cos (9) s (-g+M - - (1 - 2 cos (9)) s 

= -^(2 cos (61) sin ((f + 1)^)- sin (f0)) 

_ f sin((f+2)0) \ 

- s \ dn(fl) ; • 

and this completes the proof. 
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